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ABSTRACT 



Wind speed data is routinely acquired at the FLEET NUMERICAL 
OCEANOGRAPHIC CENTER (FNOC) in near real time from the SPECIAL SENSOR 
MICROWAVE/IMAGER (SSM/I), on board the DEFENSE METEOROLOGICAL 
SATELLITE PROGRAM (DMSP) satellite and the RADAR ALTIMETER on board the 
GEOS AT satellite. A correlation study is made of the SSM/I and ALTIMETER wind 
speed data using five data sets from areas where the DMSP satellite had crossed within 
1.1 hours of the GEOS AT satellite. Regions of ocean where the atmosphere was free 
of cloud water and precipitation have a high correlation of approximately 0.995, while 
areas with precipitation and cloud w'ater have very poor correlation. Areas of high 
water vapor content do not appear to adversely affect the high correlation between the 
SSM/I and ALTIMETER measured wind speeds. 
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I. INTRODUCTION 



Understanding the environment has been a problem mankind has dealt with for 
ages. Our knowledge of the earth was initially limited to that which we could see 
with our own eyes. Gradually our ability to sense our environment grew as new 
instruments were developed. When space satellites came onto the scene, a rich 
resource of information concerning the world around us suddenly opened. From visual 
and infrared images of the earth, the truly dynamic nature of the planet we live on 
could finally be appreciated. 

The years following the first attempts at retrieving environmental data from space 
have shown advances that parallel those made in computer technology. Devices have 
been implemented which have made use of virtually all of the electromagnetic (EM) 
spectrum. Each remote sensing device was designed for the purpose of detecting some 
form of radiation either emitted or reflected by the earth. To observe particular 
features of the atmosphere or ocean, methods were developed which would utilize 
particular frequencies to penetrate through obstacles in the path of this radiation and 
arrive at its target parameter. 

Wind speed over the ocean has proven to be a particularly difficult parameter 
to measure, but most worthwhile to pursue. Vast areas of ocean have been "data 
deserts" for years. Remote sensing of these regions using microwave frequencies can 
supply data to improve "data starved" numerical input for environmental models. 
This, in turn, should offer a better understanding of synoptic weather systems, 
improving on the visible and infrared technologies currently used today. Disseminating 
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information concerning the global surface wind field to all vessels, would enhance 
preparations for oceanic transit. Ships already underway would be better prepared for 
heavy winds and seas. 

In the case of wind speed over the ocean, two separate techniques have been 
developed. The first is the active radar device known as the GEOSAT altimeter which 
utilizes pulses of EM energy at 13.7 GHz to analyze the surface of the ocean. Though 
its prime objective was to map the earth’s geoid, its signal can also be used to measure 
surface wave height and wind speed over the ocean. The second device known as the 
SPECIAL SENSOR MICROWAVE IMAGER (SSM/I) utilizes a totally different 
approach to the problem of wind speed measurement. Instead of transmitting its own 
radiation to the earth and then receiving a returned signal for analysis, it passively 
measures upwelling radiation from the earth’s atmosphere and oceans. 

How these two methodologies compare, is the subject of this thesis. The data 
utilized in this comparison represents several areas of the world where both satellites 
were co-located within one hour of each other’s orbit. In the pages to follow, details 
of how each system operates and the physics behind the analysis of the data, w r ill be 
explored. Correlation studies made for each region will disclose similarities and 
differences in the wind speed measurements. A particularly interesting case will be 
presented where several environmental parameters have interfered with the sensors 
ability to accurately access the surface wind speed. 

Chapter two details the basic physics applicable to the SSM/I and GEOSAT 
remote sensors, as well as the operation of both systems. Chapter three will discuss 
ground truth verification studies made for each system. Here, individual strengths and 
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weaknesses are taken into account prior to the comparison study. The experiment is 
presented in chapter four, followed by conclusions and recommendations in chapter 
five. 
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II. FUNDAMENTALS OF MICROWAVE REMOTE SENSING 



The basic physics behind the remote sensing techniques examined in this thesis, 
will be discussed in this chapter. The first section will cover the propagation 
characteristics of EM radiation. This will be followed by an overview of the index 
of refraction, blackbody radiation, emissivity, and the radiative-transfer equation. This 
should provide an ample background for understanding the operation of both the SSM/I 
and the GEOSAT altimeter. 

A. PROPERTIES OF EM WAVE PROPAGATION AND ATTENUATION 

Microwave remote sensing depends on the modification of EM radiation by the 
environment. Changes take place in the phase and intensity of an EM wave as it 
passes through the atmosphere and interacts with the ocean surface. These alterations 
of the wave’s character are interpreted by the remote sensing device as variations in 
the state of the environment. We will first describe the EM wave equation and then 
develop the relationships that govern the changes that take place in the wave as it 
passes through matter. 

EM radiation consists of alternating electric and magnetic fields [Ref. 1], which 
are classified by their frequency of fluctuation. A single EM wave vector, E(x.y,z), 
can be broken down into its individual components, E„, E y , and E 2 for a simplified 
analysis (Fig.l). The propagation characteristics of a single component of an EM wave 
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Figure 1 Electromagnetic wave vector and components. 



can be described by the following equations: 



E = E 

i ° 




(1) 


k 2 = co 2 p£ - icop 


(2) 


(0 = 


2trf 


(3) 


H = 


HA 


(4) 


= e 0 e, = 


e 0 (e' - ie") 


(5) 



where: 

E y = component of the electric Field intensity 
E„ = initial electric field magnitude 
k = complex wave number 
to = frequency in radians per second 
f = frequency in hertz 
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p o = permeability of free space (4k x 10' 7 Hz/m) 
p, = relative permeability of a material (=1) 

£ = permittivity of a medium 

e Q = permittivity of free space (8.85 x 10' 12 Farads/m) 
e',e" = real and imarginary components of the dielectric constant 
For most cases in remote sensing applications, the value of p r = 1, so that 
equation (2) can be written as: 

k 2 = co 2 p 0 e 0 [e' - ia/(e c to)] (6) 

where: 

£ r = [e' - ia/(£ o 0))] the complex dielectric constant 
O = conductivity of the medium 

In free space, i.e., a vacuum, £ r = 1, a = 0, and the EM wave velocity is given 
as: 

c 0 J = l/(p„£ 0 ) (7) 

For EM w'aves traveling through matter, such as the ocean or atmosphere, the velocity 

is determined through the equation: 

c = to/k = cJ-Je, (8) 

From (7) and (8), we can define the index of refraction n, which is related to E r by 
the formula: 

n = cjc = (9) 

Equation (9) is important for understanding how the EM wave is effected as it passes 
through the atmosphere and interacts with the ocean below. Since n is related to the 
complex dielectric constant £„ radiation passing through matter will experience a 
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modification in speed and intensity. The wave will be slowed due to the resistance 
encountered by an induced electric field produced in this medium and attenuated 
through the imaginary component of n. To understand how the imaginary component 
can attenuate the wave, n must be examined in more detail. The index of refraction 
written in its component form: 

n = a + i(3 (10) 

where a and P represent the real and imaginary parts of n and are written as: 

a 2 = ( I e, I + e')/2 (11) 

P 2 = (kl - e')/2 (12) 

k I = [e' 2 + o7(e.©)T (13) 



Returning to equation (8) for the wave number k and substituting the component form 
of n into that expression, we can rewrite k as: 

k = ooa/c + iooP/c (14) 

By replacing the expression for k in equation (14) into the original expression for the 
propagating wave, equation (1), we have: 

£ _ £ e (-i(oxip/c -rot)] e Mx/c ) ( 15 ) 

= E e l<H * wt) J 

where <j) = -oo/c is known as the attenuation coefficient. With a negative value in the 
exponential, this component acts to reduce the value of the electric field as the EM 
wave passes through matter. The distance the radiation travels through some medium, 
in which the intensity of the wave is reduced by a factor of exp(-l), is known as the 
skin depth d and is given by : 

d = c/(ooP) or d = c/(2rcfp) (16) 
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For a frequency of 5 GHz, the value of d is 0.5 centimeters. The SSM/I operates 
between 19 and 37 GHz and the GEOSAT altimeter operates at 13.5 GHz, and so the 
skin depths for the instruments in this study are on the order of millimeters. 

B. POLARIZATION 

In addition to the variation of the intensity and frequency of the wave, microwave 
radiation is also characterized by its polarization. Polarization is defined as the 
orientation of the electric field vector relative to the surface it is reflected or 
emitted [Ref. 2]. Polarized radiation is classified as either horizontal, in which case 
the field is parallel to the surface, or vertical, whereby the field is perpendicular with 
respect to the surface. A change in surface conditions and viewing angle, will affect 
the vertically and horizontally polarized EM waves in different ways. The amount of 
radiation emitted from a surface (emissivity) depends on how the radiation is polarized 
(Fig. 2). The SSM/I s sensors can detect both horizontally and vertically polarized 
radiation. Certain combinations of polarized frequencies are used to allow particular 
environmental parameters to be measured separately with a minimal of interference 
from the other parameters. 

Assuming a smooth (specular) surface, the reflection of the spectral radiance, 
R, from such a surface is written in equation fonn as: 

R(X;0 r ,<}> r -7c) = L.aieAVL^e,^) (17) 

where: 

X = wavelength of the radiation 
0,4> = spherical coordinate angles 
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Incidence Angle (degrees) 

Figure 2 Polarized microwave radiation at 3 cm wavelength as a 
function of emissivity and incidence angle. 



L„L, = reflected and incident radiance respectively 



The spectral reflectance of a plane dielectric surface is determined through the Fresnel 
reflection coefficients as follows [Ref. 3]: 

R„ = [(p - cos9) 2 + q 2 ]/[(p + cos0) 2 + q 21 (18) 

R v = [(e'cosG - p) 2 + (e"cos0 + q) 2 ]/[(£'cos9+ p) 2 + (e"cos0 + q) 2 ] (19) 

where R v and R H , are the vertically and horizontally polarized components of the EM 
wave, and £ = e' - i£"; the complex dielectric constant of the surface. The quantities 
p and q are: 

p = (1/V2 )([(e' - sin 2 0) 2 +E" 2 ]" + [£' - sin 2 0]p (20) 
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( 20 ) 



p = (1/V2){[(e' - sin 2 ©) 2 +e"T + [e' - sin 2 0]p 
q = (1/V2)([(e' - sin 2 0) 2 +e" 2 ]“ - [e' - sin 2 ©]} 1 ' (21) 

where 0 is the incidence angle. 

In the case when 0=0, such as with the viewing angle of the GEOSAT 
altimeter, R H and R v are reduced to the form: 

R h = Rv = [(P - U 2 + q 2 /[(p + l) 2 + q 2 ] (22) 

= [(a - l ) 2 + p 2 /[(a + l ) 2 + P 2 ] 

where the complex index of refraction is given by equation (10). Since n 2 = e' - ie", 
the Fresnel reflection coefficient at normal incidence can take the form: 



R v = I (n - l)/(n + 1)1 2 (23) 

The spectral emissivity of a surface is related to the spectral reflectivity through 
Kirchoff’s law (20). We can relate the spectral reflectance R to the emissivity in 
both horizontally and vertically polarized forms: 



e v (0) = 1 - R V (B) 


(24) 


e H (0) = 1 - R h (0) 


(25) 



C. BRIGHTNESS TEMPERATURE 

The purpose of this section is to develop the relationship between the radiation 
emitted by the ocean’s surface and the upwelling temperature detected by the SSM/I 
known as the brightness temperature, T b . All objects, except those existing at absolute 
zero, emit radiation. An object is classified as a perfect emitter, known as a black 
body, if it radiates energy at the maximum rate predicted by Planck’s radiation law. 
As the temperature of an object increases, the total amount of energy emitted also 
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increases. Most objects, however, emit at a rate less than that of a blackbody. The 
term "grey body", is used to describe an object which emits radiation at a reduced rate. 
Variations in the temperature of a grey body relative to the assumed blackbody 
temperature of the earth can be identified and attributed to changes in the environment, 
such as wind speed, water vapor, and precipitation. 

Terms commonly used in discussing brightness temperature, or radiation 
measurement in general, are radiant energy, radiant flux, irradiance, and emittance. 
The term radiant energy, Q, refers to the amount of energy carried by the radiation 
(Joules). The radiant flux, O = dQ/dt, is a measure of rate of energy transport 
(Joules/sec). This flux defines the power available to produce a response in a 
measuring device. The irradiance, E = d<!>/dm 5 , is used in reference to the radiance 
incident on a surface. The radiance emitted from a surface is known as the emittance, 
M = dd>/dt. The character of the material from which radiation is either emitted or 
reflected is known as its emissivity. Each surface has a particular value, for example, 
the emissivity of salt water is 0.3, while that of dry soil is 0.92. It is this difference 
in emissivity that allows the SSM/I to distinguish between land and sea surfaces. 

The GEOSAT relies on the degree of reflectivity to detect changes in the ocean’s 
surface. The GEOSAT and SSM/I observations are linked by Kirchoff’s law, which 
relates emissivity and reflectivity in the following manner: 

e = 1 - r (26) 

where r is the reflectivity of the surface and e is the emissivity. 

The ratio of spectral brightness emitted in a direction (0,(j)) by a real body 
relative to that which is radiated by a blackbody is the spectral directional emissivity 
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e of that body, where £ is a function of wavelength A. and spherical angles <j> and 0. 
The microwave emission from the ocean surface is given by eB v , where B v follows 
the Rayleigh-Jeans approximation of Planck’s Law. Planck’s Law states that the spectral 
brightness per unit of frequency is expressed as: 

B v = (2hv 2 /c 2 ) (l/[e ,hvAkT) ' - 1]) (27) 

where: 

k = 1.38 x 10- 23 J/K, Boltzman’s constant 
h = 6.63 x 10 -34 J s, Planck’s constant 
T = temperature in degrees Kelvin 
v = frequency of radiation in hertz 
c = 3 x 10* m/s, speed of light in a vacuum 
The Rayleigh-Jeans approximation is made at microwave frequencies, where hv/(kT) 
« 1, so that: 

e lhv/<vr>i = ! + hv / (kT ) 

Substituting the above relationship into equation [18] gives 

B v = (2v 2 k/c 2 )T (28) 

Equation [19] states that the spectral brightness is proportional to the temperature of 
the body emitting radiation. The radiance received by a microwave sensor is known 
as the brightness temperature. The brightness temperature T b , is given by eT and will 
be discussed in a later chapter. 

The following section will cover the topics of absorption and the interaction of 
radiation with surfaces. These factors contribute to the variation in brightness 
temperatures as measured from satellites. 
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I). 



ABSORPTION AND St'ATTi RING OF RADIATION 



Radiation emitted by the ocean surface must pass through the atmosphere once 
before being detected by the SSM/I, and twice for the GEOSAT altimeter. 
Consequently, interactions which take place between the atmosphere and the radiation 
must be understood in order to interpret a retrieved signal. In the atmosphere, energy 
is both scattered and absorbed. For microwave frequencies, the primary absorber is 
water vapor, while scattering is predominantly due to rain. Any instrument used to 
measure wind speed over the ocean must utilize frequencies where the effects of 
atmospheric absorption and scattering are small. The GEOSAT altimeter frequency was 
chosen such that atmospheric effects would be small. In addition to observing the 
ocean surface, the SSM/I measures atmospheric parameters such as water vapor and 
precipitation. This allows the SSM/I to adjust for atmospheric parameters which would 
interfere with the measurement of wind speed. Questionable wind speeds, based on the 
presence of these interfering parameters, are identified at ground stations when the data 
is analyzed. 

E. THE RELATION OF SURFACE ROUGHNESS TO EMITTED RADIATION 

The roughness of the ocean causes the reflected and emitted EM waves to be 
sent out in numerous directions. Remote sensors of microwave radiation must be able 
to account for this variation in directionality, hence the viewing angle becomes an 
important factor. The received radiation is also a function of the dielectric constant 
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of the surface. These properties can he related if we, (1) assume that the surface is 
made up of many small facets, (2) examine only a single facet, or plane dielectric 
surface and apply the Fresnel Reflection Coefficients to that surface to obtain the 
spectral reflectance of that surface. Finally, (3) integrate over all the facets in the area 
observed by the sensors to obtain the spectral reflectance for that surface. 

The solution to the Fresnel equations requires that either the complex dielectric 
constant, e t , or the index of refraction of the surface as a function of frequency be 
known. In the case of the ocean, the dielectric constant is a well known parameter 
through the work of Lane and Saxon [Ref. 4], The dielectric constant is represented 
by the Debye equation: 

£, = £„ + [(£,£_)/( 1 + (ion) 1 - 0 - ia/coej (29) 



where: 

co = 2jtf, radiation frequency 

£„ = dielectric constant at infinite frequency 

£, = the static dielectric constant 

t = the relaxation time 

a = the ionic conductivity 

a = an empirical constant 

The parameters £„ t, and c are all functions of the material which is reflecting or 
emitting the radiation to be received by microwave sensors. The conductivity a of 
the sea water increases the complex portion of the dielectric constant, however, it is 
important to note that the influence of sea salts on the emissivity of the ocean is small 
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for frequencies greater than 5 GHZ. Since both the SSM/I and GEOSAT sensors 
operate above 5 GHZ, the salinity of the sea water has no effect. 

The above description for spectral reflectance and emittance from a single facet 
can now be applied to the entire rough surface. This is accomplished by integrating 
over all the facets and taking a statistical average. 

F. THE RADIATIVE TRANSFER EQUATION 

The radiative transfer equation contains terms which relate radiant energy from 
the surface and atmosphere to the radiant energy observed by a satellite. Before 
energy can reach a satellite sensor from the ocean surface, the upwelling radiation is 
attenuated and joined by additional emission from the atmosphere, so that the observed 
temperature detected by the SSM/I is given by: 

T = T a (1 - t) + t[eT. +(1 - e)(l - t)TJ (30a) 

or 

T = et[T, - T a ( 1 -t>] + T A (1 - t 2 ) (30b) 

where, T A is the effective temperature of the atmosphere, t = e' '\ x is the optical 
depth, T, is the sea surface temperature, and e is the emissivity. The brightness 
temperature of the radiation incident on the ocean surface by the atmosphere is 
T A (l-t). The amount that is reflected by the sea surface is written as (1 - e)T A (l - 
t), while the amount emitted is eT,. It follows that the total upwelling brightness 
temperature from the sea surface is eT, + (1 - e)(l - t)T A . The signal from the 
surface is eT,. The other terms are considered noise with regard to surface wind 
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speed analysis, but are useful when detecting water vapor, liquid water, and 
precipitation. If the optical depth x = 0, then the atmosphere is transparent and allows 
just the brightness temperature of the ocean surface to be seen by the sensor. Equation 
(30b) would then be reduced to simply T = eT,. This the ideal case for measuring wind 
speed over the ocean. If, on the other hand, x » 0, then T = T A , i.e., the SSM/I 
would be looking at the environmental changes in the atmosphere only. The 
frequencies and the method used by the SSMI and GEOSAT to minimize atmospheric 
interference will be discussed in sections devoted to instrument operation. 

G. THE TRANSFER OF MOMENTUM FROM THE ATMOSPHERIC SURFACE 
LAYER TO THE SEA SURFACE 

Waves and foam generated by the wind contribute to the roughness of the sea 
surface and effect both microwave emission as measured by the SSM/I and microwave 
scattering as measured by the GEOSAT altimeter. The roughness has the effect of 
altering the incidence angle of radiation being emitted or reflected and changing the 
emissivity of the sea surface itself. The roughness of the sea surface is an indication 
of the transfer of momentum into the sea from the turbulent atmospheric surface layer. 
An extrapolation process is performed to relate the wind speed in the atmospheric 
surface layer to the sea surface roughness. This section will briefly describe the 
transfer of momentum from the atmospheric surface layer to the sea surface. 

Businger [Ref. 5] describes the formulation of the turbulent transfer of momentum 
from the atmospheric surface layer to the sea surface as a semi-empirical relationship. 
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The appropriate equation can be developed by utilizing the turbulent transfer coefficient 
K in the equation for momentum transfer, where: 

F m = -pK9U/a Z (31) 

Here, p is the density and U is the mean horizontal wind speed. Businger defines 
the atmospheric surface layer as the lower 20 meters of the atmosphere. Here the 
interaction with the surface is strong and the layer adjusts to the surface conditions 
rapidly. We can therefore assume a quasi-steady-state. The momentum transferred 
through to the ocean surface exerts a drag force per unit area on that surface known 
as the surface stress t. By Newton’s third law, the momentum equation is given as: 

F m = pu'w' = -T (32) 

The friction velocity u. is commonly used as a scaling velocity, where: 

u. = (33) 

Equation (31) can now be written as: 

K.a’U/dz = u. 2 (34) 

where u. is independent of height. To generate a wind profile from the surface to a 
height of 20 meters, an assumption must be made about K^. By dimensional analysis, 
has the dimensions of velocity times length. We use u. to scale the velocity and 
z to scale the length. Therefore, we assume that: 

= ku.z (35) 

where, k is a constant of proportionality, the von Karmann constant. Combining 
equations (34) and (35), we have: 

9U/dz = u./(kz) (36) 
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To avoid an infinite shear at the surface, when z = 0, a finite surface roughness z D is 
introduced into the equation. This now becomes: 

d\J/dz = u./(k[z + zj) (37) 

Upon integration we have the form: 

U/u. = 1/k ln[(z + z 0 )/z 0 ] (38) 

which is the logarithmic profile for the wind from the sea surface through the 
atmospheric surface in a neutrally stable atmosphere. This is the common profile used 
to extrapolate surface wind speed to a height of 10 or 20 meters, which are the typical 
heights at which sea buoys and ships measure wind speed at sea. The neutral case is 
used primarily for its simplicity and also for its adequate approximation of the surface 
layer wind field. Figure (3) [Ref. 6] represents a semi-log plot of wind speed from 
the surface through the atmospheric boundary layer for neutral, stable and unstable 
conditions. It would be impractical to measure the stability of the atmosphere for all 
wind speed measurements. Fortunately, the neutral profile of wind speed approximates 
the atmosphere in most situations. 

Fig. 3 is derived from setting equation (37) equal to a wind shear function, <t> M , 



such that: 


0> M = (k[z + 


zj/u.) dU/dz 




where: 










<D M = 1 + [ 4.7 z/L ], 


for z/L > 0 (stable) 


(39a) 




= 1 , 


for z/L = 0 (neutral) 


(39b) 




<D m = [1 + [ 15 z/L ]r. 


for z/L < 0 (unstable) 


(39c) 



and L is the Obukhov length. 
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Height 




Figure 3 Typical wind speed profiles versus static 
stability in the surface layer. 
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III. MEASURING WIND SPEED OVER THE OCEAN WITH THE SSM/I 



The SSM/I is a microwave radiometer which measures upwelling radiation from 
the earth’s surface at five frequencies and two polarizations. Using frequencies 
associated with atmospheric windows, the SSM/I measures upwelling thermal 
microwave radiation from the ocean surface. This radiation is called the brightness 
temperature. From equation (30b), where T b = eT, for an atmosphere with 100% 
transmittance, the brightness temperature is a function of the emissivity, e, and 
temperature of the sea surface, T,. The sea surface temperature is relatively constant 
compared to the change in emissivity. An increase in T b is thus caused by an increase 
in the sea surface emissivity. A change in surface roughness and thus wind speed, will 
change the emissivity of the surface. Thus, a particular wind speed at the sea surface 
can be related to a certain brightness temperature. 

The SSM/I utilizes the 19.35H, 22.235V, 37H, and 37V GHz channels for 
measuring wind speed, where H and V represent horizontal and vertical polarizations 
respectively. The 19.35H, 37H, and 37V GHz channels are "windows" in which there 
is nearly 100 % transmission of upwelling radiation through the atmosphere. The 
22.235 GHz channel is in the center of a water vapor absorption line and is used to 
measure water vapor concentration. Detection of water vapor is vital to the 
determination of the validity of the wind speed detected, since the 19.35H, 37V, and 
37H channels are subject to some small amount of atmospheric interference which 
require that an atmospheric correction be applied. 
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The algorithm used by the SSM/I for measuring wind speed assumes the 
toughness of the ocean surface can be treated as a conglomeration of many tiny 
emitting planes whose angle of inclination is altered by the force of the wind. The 
variation in the sea surface slope is used in calculating the emissivity of the surface 
by assuming a Gaussian slope distribution of the angles of the slope. A change in 
the distribution of these angles corresponds to a change in the root mean square slope 
of the sea surface and distribution of the emissivities associated with the variance of 
sea surface slopes. The resulting change in emissivity causes a change in brightness 
temperature through the Fresnel equations (Eqns.22 and 23). The equation used to 
describe the dependence of derived wind speed on the surface slope variance, o J 
[Ref. 7] , is as follows: 

cr(f) = (0.3 + 0.02f) (0.003 + 0.48U 2O ) for f < 35 GHz (40a) 

c J (f) = 0.003 + 0.48Ujo for f > 35 GHz (40b) 

Here, f is defined as the frequency in GHz and U 20 as the wind speed measured in 
meters/second at a 20 meter height. 

Foam generated by breaking waves occurs at wind speeds greater than 7 m/s. 
The affect of foam on the microwave emission from the sea surface is not well 
understood and there are varying theories in literature. Swift [Ref. 8] cites an 
experiment over the Salton sea using a 19 GHz radiometer to measure wind speed 

over a foam covered sea. Swift states that the results were tens of degrees higher than 

theory and could be attributed to the presence of foam. Wilheit [Ref. 9) treated foam 
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as a relatively non-reflective surface and eliminated it from the radiative transfer 
equation. The reflectivity that does occur comes from the dielectric interface around 
the foam bubbles. Hollinger et al [Ref. 10] treat foam as "a substance which partially 
obscures the ocean surface in a manner independent of polarization but dependent upon 
frequency". He defines the foam fraction, F f , as 

F f = 0.006 (U„ - 7) (1 - c (ns ) for U M > 7 m/sec (41a) 

F f = 0 for Ujp < 7 m/sec (41b) 

This equation is developed from the observation that the sea surface brightness 
temperature increases linearly for wind speeds in excess of 7 m/sec, and that there is 

an absence of foam development for speeds less than or equal to 7 m/s. The foam 

fraction is also used in the empirically derived SSM/I wind speed algorithm to account 
for the change in emissivity that occurs with the presence of foam. 
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IV. MEASURING WIND SPEED OVER THE OCEAN WITH THE GEOSAT ALTIMETER 



The primary mission of the GEOSAT satellite was not wind measurement, but 
measurement of the geoid and sea surface height. However, wind is computed as a 
by-product and with an accuracy of +/- 2 m/s [Ref. 10]. To measure wind speed, the 
GEOSAT satellite utilizes an active radar operating at a frequency of 13.9 GHz called 
an altimeter. This instrument transmits a pulse of microwave energy to the ocean 
surface and measures the return signal. The altimeter [Ref. 11] uses an automatic gain 
control (AGC) loop to normalize the amplitude of the return signal from the sea 
surface. The AGC is a measure of the backscatter cross section of the sea surface, a , 
and is related to the wind speed. The greater the wind speed, the lower the cross 
section, thus the weaker the return signal from the ocean surface. The backscatter 
cross section is obtained from the following equations [Ref. 12]: 

AGC [dBm] = 10 In (PJPj (43) 

where: 

P r = the received power 

P 0 = the reference power (1 mW) 

The received power P r , is defined in terms of antenna gain G, transmitted power P„ 
the operating wavelength X, the distance to the sea surface R, and the scattering cross 
section a . The relation appears as follows: 

P r = [AV(4 tx>] P, G a / R 4 dA (44) 
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Assuming that the sea surface slope is both gaussian and isotropic in its probability 
distribution, we can define a as follows: 

a = [ I R(0) 1 2 / S 2 ] sec 4 © exp [ -tan 2 © / S 2 ] (45) 

where: 

0 = angle of incidence 

R(0) = Fresnel reflection coefficient for the air-sea interface at normal incidence 

The Cox and Munk [Ref. 13] relationship for average upwind/crosswind mean 
square sea surface slope a 2 and wind speed is given as: 

a 2 = a + 6 U (46) 

where, U is the wind speed in m/s and a, B are empirical constants of proportionality 
which depend on the probability density of the sea surface slope. By making the 
assumption of an isotropic Gaussian surface, the separate effects of the upwind 
component and the crosswind component are ignored [Ref. 14], By ignoring wind 
direction when deriving a relationship between wind speed and surface slope, error in 
the surface wind speed will be introduced. 

If we let 0 = 0 , since we are observing at the nadir angle, the equation for a 
becomes: 

a = I R(0) 1 2 /(a + BU) (47) 

The equation for the received power is 

P r = a 2 /(47t) 3 ) P, (G 2 /h 4 ) a A (48) 

where: 

h = satellite altitude 
A = area illuminated by the altimeter 

This relationship is used in the Brown algorithm for the actual calculation of wind 
speed [Ref. 15] at the Fleet Numerical Oceanographic Center. Brown’s algorithm is 
a result of a comparison of wind speed data received from sea buoys to normalized 
backscattered cross sections derived from altimeter power measured on board the 
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GEOS-3 satellite and applied to GEOSAT. The algorithm is presented in two stages. 
The first section is as follows: 



W, = exp [(io- <02,Io710 > - B)/A] 



(49) 



where: 

A = 0.080074 B = -0.124651 for a* < 10.12 dB 

A = 0.039893 B = -0.031996 for 10.12 < a < 10.90 dB 

A = 0.015950 B = 0.017215 for a" > 10.9 dB 

In the above equation, the backscattered cross section is given as decibels relative to 

1 m 2 . 

A second equation is used to complete the derivation of surface wind speed. This 
was necessary to adjust for the skewness in the distribution of the differences between 
the buoy and altimeter wind speed measurements. The second part of the algorithm is 
as follows: 



5 n 



2 A, W, 

n*l 


, for W, 


< 16 meters per second 


(50a) 


w„ 


for W, 


> 16 meters per second 


(50b) 



where: 

a, = 2.08779900000 
a 2 = -0.36499280000 
a 3 = 0.04062421000 
a 4 = -0.00190495200 
a, = 0.00003288189 

U I0 is the w'ind speed at a height of 10 meters. This algorithm appears to measure 
the wind to within +/- 2 m/s. 
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V. GROUND TRUTH STUDIES 



A. GEOS AT Wind Speed Ground Truth Study 

Shuhy [Ref. 16] made a ground truth study for the wind speeds obtained by 

GEOSAT using National Data Buoy Center buoy data. To avoid possible interference 

from land, most of the buoys used were located in deep water away from any 

coastlines. Shuhy chose to restrict the comparison to satellite passes within 80 

kilometers and 90 minutes of a ground truth point. The data was classified as either 

excellent, good, or fair based on the following criteria: 



CLASSIFICATION 


DISTANCE FROM 
GROUND TRUTH 


TIME DIFFERENCE 
IN OBSERVATIONS 


excellent 


25 KM 


30 min 


good 


25 - 60 KM 


60 min 


fair 


60 - 80 KM 


90 min 



No data outside 80 kilometers or 90 minutes from crossing were used. The results of 
the comparison are seen in Fig. 4. Shuhy concluded from regression analysis that the 
sub-categorization described above yielded no significant differences. A total of 256 
data points were used with an rms error of 1 .54 m/s. The comparison demonstrated that 
GEOSAT wind speeds are higher than the ground truth for speeds less than 6 m/s and 
for speeds greater than 6 m/s, GEOSAT winds are lower than those obtained on the 
surface. Shuhy notes, however, that there is significant scatter in the data and that 
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Figure 4 Geosat vs ground truth winds at a height of 10 meters. 



spatial variability must be examined to determine which comparisons might be invalid 
due to high variability between buoy and satellite wind speeds along the satellite track. 

When ground truth data is compared to off-nadir GEOSAT data, there appears to 
be a correlation between the magnitude of the angle and how well the GEOSAT and 
buoy wind speed compared. For angles greater than 1.1 degrees, the correlation was 
poor. Angles less than 0.4 degrees off nadir exhibited good correlation in the 
GEOSAT/buoy wind speed comparisons. In this study, the nadir angle was not 



27 



DMATRIX WINDS (M/S) 



available as the GEOSAT data used was preprocessed by Naval Ocean Research and 
Development Activity (NORDA), however all data with off nadir angles greater than 
1' is not transmitted to FNOC. 

B. SSM/I WIND SPEED GROUND TRUTH STUDY 

The SSM/I algorithms have undergone four revisions since the system was put 
into space in late 1987. The study completed by Hollinger and Poe [Ref. 18] includes 
600 SSM/I and buoy data sets. Comparisons of wind speed measurements were made 
within 30 minutes and 25 kilometers of each other (Fig. 5). 




BUOY WINDS (M/S) 

Figure 5 SSMI wind speed comparison with ground truth from moored 
buoys (Hollinger and Poe). 
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The SSM/I winds are extracted from a data set by the "D-matrix" method. This 
is an empirical/statistical approach that chooses the most probable atmospheric and 
surface conditions that would produce a characteristic set of measured brightness 
temperatures. The "D-matrix" contains the reference brightness temperatures for the 
SSM/I. The SSM/I retrieved brightness temperatures are compared to the "D-matrix" 
and from this comparison the ocean and atmospheric conditions that were observed 
by the SSM/I are determined. The overall scatter of data points represents a +\- 2 
m/s standard deviation between 3 and 15 m/s. At this time, there is insufficient data 
for wind speeds greater than 15 m/s to make an accurate comparison of SSM/I and 
buoy wind speeds. 

Interfering factors that could reduce the accuracy of the analysis are heavy 
concentrations of water vapor and the presence of cloud water and rain. In this study, 
areas exhibiting these conditions were included in the data to observe their influence 
on retrieved wind speeds. 

Positioning error was discovered by as much as 20 to 30 km for shift in pixel 
latitude and longitude coordinates. However, Hollinger and Poe state that this did not 
interfere with the overall validation of the wind speed algorithm. 
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VI. SSM/I AND GEOSAT ALTIMETER WIND SPEED COMPARISON 



A. DATA ACQUISITION 

The data base used was acquired from Fleet Numerical Oceanography Center 
(FNOC) in Monterey, California. All data obtained for this study was gathered in 
near-real time using a personal computer (PC) with a 1200 baud modem. Co-located 
passes of the DSMP and GEOSAT satellites within one hour of each other were 
required. The position and time of ascending node for the orbits of each satellite 
were calculated well in advance. These orbit tables were updated every two weeks to 
account for errors in calculating the exact times of ascending node crossings. 
Adjustments had to be made to the GEOSAT orbit table to account for the periodic 
maneuvering of the satellite that would keep it in a 17 day repeat orbit pattern. To 
avoid any microwave interference from land stations, areas over the open sea were 
chosen. 

Several problems were encountered in obtaining data in this manner. The average 
length of time to download a single SSM/I data set was at 1200 baud was 
approximately 0.5 to 1 hour, depending on its size. Near real time SSM/I and 
GEOSAT data at FNOC is accessible for 6 hours or less before it is removed from 
memory. As a result, many data sets obtained were incomplete as the data had been 
removed during the downloading process. I was unable to determine at what times a 
particular orbit would be resident at FNOC and for what length of time it would be 
resident. I consequently worked with the data available at the time. Another difficulty 
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arose in the form of the GEOSAT data format. The GEOSAT data obtained from 



FNOC was pre-processed by NORDA in Bay St. Louis, Mississippi. The format of 
the GEOSAT data base included only the latitude, longitude, wind speed, and 
significant wave height. Raw data would have been ideal, so that various wind speed 
algorithms could have been tested. Occasionally incomplete data sets were obtained 
where wind speeds were missing for large areas of an orbit. This was due to data 
transmission problems occurring between NORDA and FNOC. Despite these 
difficulties, five data sets were obtained from three different regions of the globe: the 
North Atlantic, the Gulf of Alaska, and a region in the North Central Pacific. 

B. DATA FORMATS AND PROCESSING 

This section describes the various formats in which the SSM/I and GEOSAT 
data are found. In addition, the process in which raw SSM/I infonnation is decoded 
and then combined with the preprocessed GEOSAT wind speeds to form the final data 
base will be discussed. 

1. SSM/I Data Base 

There are four forms of data which are on file for each orbit made by the 
DSMP satellite. All are encoded in hexadecimal format and must be decoded to 
decimal before processing. 

The Orbit Summary File (SUM) contains information on the latitude and 
longitude of each end of the sensor swath for entire satellite orbit. Included with this 
is a history of sensor health, i.e., whether or not any of the sensors on board 
experienced any unusual fluxes in radiation, including the hot calibration standard. 
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reference voltages, and the automatic gain control From this file, the exact time during 
which the DMSP orbit would coincide with a GEOSAT crossing could be determined. 
This fixed the amount of data needed to be downloaded to cover as much of the 
GEOSAT track as possible. 

Raw brightness temperatures for each sensor, are initially stored in a file known 
as the Sensor Data Record (SDR). An Antenna Pattern Correction (APC) is applied 
to the brightness temperatures to correct for spurious energy received in the side lobes. 
Corrected brightness temperatures are then stored in the Temperature Data Record 
(TDR). It is in this format that brightness temperatures can be used to derive 
geophysical data such as precipitation rates, ice edge, ice age, water vapor, and wind 
speed. The geophysical data are then stored in the Environmental Data Record (EDR). 
This is the level of data used in this study. 

2. GEOSAT Data Base 

GEOSAT altimeter data is stored on board the satellite for approximately 
12 hours and then transmitted to a receiving station at the Johns Hopkins Applied 
Physics Laboratory (APL) where it is archived on tape. Three data products are 
produced from the raw data; The Sensor Data Record, the Naval Ocean Research and 
Development Activity (NORDA) Data Record, and the Waveform Data Record. The 
first two contain information about the measured altitude, significant wave height, wind 
speed, automatic gain control, and corrections for satellite and instmment errors. The 
last record contains waveform data used primarily for ice studies. 

FNOC stores GEOSAT data received from NORDA in half hour blocks 
which contain latitude, longitude, wind speed, and sea surface height information. 
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The data is pre-processed and does not contain the raw data necessary to verify the 
algorithms used to produce the wind speeds and corrections made to the data set. 
Reliance was placed on the ground truth comparison made with the GEOSAT wind 
speed calculations to establish the validity of the wind speed data. 

3. Processing GEOSAT Altimeter and SSM/I data sets 

Several Fortran programs were used to convert the hexadecimal format of 
the SSM/I data records to decimal. A program called COMPARE was written to 
combine the measurements from both the GEOSAT and SSM/I EDR data bases for 
each co-located satellite pass. SSM/I wind speeds measured within a 0.25 degree 
radius of a GEOSAT data point were averaged together and a mean value and standard 
deviation were obtained. A new data set was then produced which contained the 
GEOSAT wind speeds and the mean values of the SSM/I wind analysis from which 
correlation studies were made. A graphical representation of the two data sets were 
produced. The SSM/I EDR also supplied the rain, cloud water, and water vapor 
analysis necessary to determine the environmental conditions that existed during the 
wind speed measurements. 

C. DATA ANALYSIS AND RESULTS 

A comparison of all the compiled data sets for the SSM/I and GEOSAT altimeter 
wind speed analysis is seen in Fig. 6. The correlation for the compiled data sets is 
0.703. From this graph, four distinct regions are seen. The region "A" extends from 
2.5 to 9 m/s. Here, the SSM/I demonstrated good correlation with the GEOSAT 
altimeter. Wind speeds between 5 and 8 m/s represent compiled data from four 
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different orbits, however, wind speed from 2.5 to 5 m/s was derived from a single data 
set, which resulted in a region of sparse data. The SSM/I is limited to a minimum 
detectible wind speed of 3 m/s. This limitation is due to the lack of sensitivity to 
changes in emissivity for wind speeds between 0 and 3 m/s. The region greater than 
8 m/s demonstrated very poor correlation between the SSM/I and the GEOSAT 
altimeter. This area is divided in to three separate regimes. The region "B" 
corresponds to the bulk of data which is greater than 9 m/s. Here, the correlation 
appears to be good, but the SSM/I data are notably biased to the higher wind speeds. 
On either side of region "B", there are outlying points. Region "C" is identified by 
the extremely high values of wind speed of the SSM/I for GEOSAT winds from 8 to 
12 m/s. Region "D" is represented by the high values of the GEOSAT wind speeds 
relative to the SSM/I values from 11 to 14 m/s. 

As both systems have successfully verified wind speeds in the regime between 
8 to 15 m/s, the regions where the two measurements differed were examined in detail. 
The water vapor, cloud water content, and rain measured by SSM/I were used to 
determine if any interference due to these substances may have caused the bad 
correlation in the high wind speed areas. This was done for each individual data sets 
to determine if regional differences may have made a contribution. Data sets are 
identified by the corresponding SSM/I orbit number, latitude and longitude of the 
region, as well as the date and time difference between the SSM/I and GEOSAT co- 
located orbits. 
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COPIIEO DATA SETS 




1 - 



w l 1 1 1 

0 2 


T 1 1 1 1 1 1 1 I 1 1 I 1 I 1 

i 6 8 10 12 H 16 18 

ALTlACTEfi WIK) StED C^S) 


Figure 6 SSMI versus GEOSAT wind speeds in meters per second. 
Line number 1 represents the actual correlation of 0.703 and line 2 
represents a reference correlation of 1.0. 

1. Orbit 9138, 27 March, 1989 


The scan area of the 


SSM/I encompassed a region bounded by the following 


coordinates: 




Latitude 


Longitude 


32.50°N 


226.00°E 


32.50°N 


242.20°E 


17.00°N 


226.00°E 


17.00°N 


242.20°E 



35 



The time difference between the GEOSAT and SSM/I passes is no greater than 0.61 
hours. The track of the GEOSAT through the SSM/I swath begins at 20°N,231.8°E 
and ends at 29.5°N,227.4°E. The wind speeds recorded in this region ranged from 
2.8 to 7.9 m/s. The correlation for this pass is 0.925 (Fig. 7). The difference in 
wind speed for the individual satellites demonstrates a maximum of 1.9 m/s (Fig. 8). 
The wind field, as observed by the SSM/I, is depicted in Fig. 9. The entire area of 
co-located orbits reveals a region of high water vapor content of approximately 60 
kg/m 2 (Fig. 10). No cloud water or precipitation is present in this region. Since the 
correlation of the wind speed for orbit 9138 is so high, one could safely deduce that 
a water vapor content as high as 60 kg/m 2 does not adversely affect the remote sensing 
of wind speeds by microwave instruments. 

2. Orbit 8934, 13 March, 1989 

The scan area of the SSM/I encompassed a region bounded by the following 
coordinates: 



Latitude 


Longitude 


57.94°N 


214.33°E 


54.44°N 


192.55°E 


47.39°N 


216.97°E 


44.57°N 


199.38°E 



The time difference between the GEOSAT and SSM/I passes is within 1.01 hours. 
The track of co-located orbits begins at 39.5°N,217.8°E and ends at 27.8°N,21 1 .5°E. 
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55M/ I VS GE05AT WIND SPEED 




Figure 7 SSM/I vs altimeter wind speed. Line number 1 represents the 
actual correlation of 0.925 and lute 2 represents a reference correlation 



of 1.0. 
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Figure 8 SSM/1 and GEOSAT wind speed vs. latitude, orbit 9138. 
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Figure 9 SSM/1 wind speed field (m/s): orbit 01 38. L.ine through field represents 
GEOSAT flight path. 
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Figure 10 SSM/I water vapor field (kg/m 2 ); orbit QI38. The line through the field 
represents the GEOSAT flight path. A high correlation of wind speeds exist 
between the GEOSAT altimeter and SSM/I despite high water vapor concentration. 
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A break in the track was due to an interruption in communications during the transfer 
of data over a modem. Data were consequently lost between 31°N and 32.5°N 
latitudes. The statistics for this orbit revealed a correlation between the GEOSAT 
altimeter and SSM/1 wind speeds of 0.616. This data set spans a range of wind speeds 
from 5 to 9.5 m/s (Fig. 11). The data points appear to be slightly biased toward the 
higher wind speeds for the GEOSAT altimeter. The plot of wind speed versus latitude 
(Fig. 12) reveals that both satellites follow the trend of increasing and decreasing wind 
speeds. The largest difference in wind speeds (approximately 2.5 m/s) appears at 
29.2°N, 32.5°N, 35.5°N, and 38°N latitudes. The images of wind speed (Fig. 13), 
water vapor (Fig. 14), and precipitation (Fig. 15) show that the track of the co-located 
orbit does not pass through any interfering atmospheric substances and that there were 
no unusual gradients in the vicinity. Overall, orbit 8934 demonstrated good correlation 
between the GEOSAT altimeter and the SSM/I. 

3. Orbit 8058(7), 10 January, 1989 

The scan area of the SSM/I encompassed a region bounded by the following 
coordinates: 



Latitude 



Longitude 



59.73°N 



217.54°E 



57.84 C N 



197.63°E 



51.88°N 



223.64°E 



48.89°N 



205.01°E 



40 



55M/I VS. ALTIMETER WIND 5PEED 
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Figure 11 SSM/I vs. GEOSAT WIND SPEED: ORBIT 
8934. Line 1 represents the linear regression line. Line 
2 is a reference slope of 1.0. 
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Figure 12 SSM/I and GEOSAT winds vs. latitude: orbit 
8934. 
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Figure 13 SSM/I wind speed field (m/s): orbit 8934. Line through the wind 
field represents the GEOSAT flight path. 
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Figure 14 SSM/I water vapor field (kg/m J ): orbit 89?4. Line through field 
represents the GEOSAT flight path. 
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Fimire 15 SSM/I precipitation field (mm/hr); orbit 8934. Line through field 
represents the GEOSAT flight path. 
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The time difference between the GEOSAT and SSM/I passes is no greater than 1.10 
hours. The track of the GEOSAT through the SSM/I swath begins at 49°N,219°E and 
altimeter for this particular orbit. By inspection, the correlation between the two 
instruments is poor. The calculated correlation is 0.208. Fig. 17 illustrates the 
variation in wind speeds going from latitude 49°N to 57.5°N. Beginning at the most 
southern point, the difference between the SSM/I and GEOSAT altimeter grow from 
approximately 1 m/s to as much as 3.5 m/s. The corresponding SSM/I wind field is 
represented in fig. 18. Examination of the SSM/I water vapor field (Fig. 19), shows 
the orbit tracks cross through a gradient of water vapor which changes from 
approximately 10 to 30 kg/m 2 over a distance of about 50 nm. The track continues 
through an area of saturated water vapor, i.e., greater than 60 kg/m 2 and then exits into 
a region of approximately 30 kg/m 2 . The crossing of the gradient at 51°N coincides 
with the rapid increase in the SSM/I wind speed, while the altimeter remains relatively 
constant. Examination of the SSM/I precipitation field (Fig. 20) reveals precipitation 
rates from 1 to 4 mm/hr at 51°N and from 53.5°N to 56°N, along the co-located orbit 
track. From 56°N to the end of the track, cloud water is present up to 0.125 kg/m 2 
(Fig. 21). This suggests that in regions of cloud water and precipitation, SSM/I winds 
are greater than GEOSAT values. 
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SSM/ I VS ALTIMETER WIND SPEED 




Figure 16 SSM/I vs. GEOSAT wind speed: orbit 80587. Line 1 
represents an actual correlation of 0.208, while line 2 represents the 
reference correlation of 1.0. 
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Figure 17 SSM/1 and GEOSAT wind speed vs. latitude: orbit 80587. 
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Figure 18 SSM/I wind field (m/s); orbit 8058(7). Line through field represents the 
GEOSAT flight path. 
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Figure 19 SSM/T water vapor field (kg/m J ): orbit 8058(7). Line through field 
represents the GEOSAT flight path. 
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Figure 20 SSM/I precipitation field (mm/hr): orhit I.ine through wind 

field represents the GEOSAT flight path. 
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Fipuro 21 SSM/I cloud water field fkg/m 3 ): orbit 8058(7). I ine through field 
represents the GEOSAT flight path. 
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4. Orbit 8058(6), 10 January, 1989 



The scan area of the SSM/I encompassed a region bounded by the following 
coordinates: 



Latitude 


Longitude 


51.88°N 


223.64°E 


48.89°N 


205.01°E 


42.14°N 


225.01°E 


39.60°N 


210.1 1°E 



The time difference between the GEOSAT and SSM/I passes is within 1 .00 hours. 
The track of the GEOSAT through the SSM/I swath begins at 40.8°N,224.6°E and 
ends at 49.5°N,218.6°E. The corresponding graph of SSM/I vs GEOSAT wind fields 
is presented in Fig. 22. This data set spans a range of wind speeds from 7 to 12 m/s. 
The correlation in this graph is good below 8 m/s and then the SSM/I appears to read 
slightly higher then the GEOSAT altimeter. Overall, the correlation is calculated to 
be 0.935. A graph of wind speed versus latitude (Fig. 23) shows that both instruments 
w'ere w'ithin 1 m/s from 41°N to about 46.5°N and afterwards the difference began to 
open to about 2 m/s near 48.5°N. Examination of the SSM/I wind field (fig. 24) 
shows that the location where the difference in wind speed begins to grow occurs 
around the same area that the wind speeds begin to increase rapidly. The SSM/I water 
vapor field (Fig. 25) does not demonstrate any significant gradient which would 
interfere with the measurement of wind speed or any area on increased water vapor. 
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SSM/ I VS ALTIMETER WIND SPEED 



0R9IT 80586 




GEOSAT WIND SPE60 C M/ S} 

Figure 22 SSM/I vs. GEOSAT wind speed: orbit 80586. Line 1 
represents an actual correlation of 0.935, while line 2 represents a 
reference correlation of 1.0. 
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Figure 23 SSM/I vs. GEOSAT winds vs. latitude: orbit 8058(6). 
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Figure 24 SSM/1 wind speed field (m/s); orbit R058(b). I ine represents the 
GEOSAT flight path. 



53 




Fii>nrc“ 25 SSM/1 water vapor field (kg/m ? ): orbit 805R(6h I ine through field 
represents the GEOSAT flight path. 
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Figure 26 SSM/1 precipitation field (mm/hr): orbit 8058(6). Tlie line through the 
field represents the GEOSAT flight path. 
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There is, however, a gradual increase in water vapor content in the northeast quadrant 
of the field from 15 to 30 kg/m 2 . The precipitation field (Fig. 26) shows precipitation 
to the northeast of the GEOSAT-SSM/I track, but again, there is no evidence of 
precipitation occurring along the track. Overall, there appeared to be no interfering 
atmospheric substances in the region of the co-located orbits. 

5. Orbit 7925, 2 January, 1989 

The scan area of the SSM/I encompassed a region bounded by the following 
coordinates: 



Latitude 


Longitude 


79.98°N 


335.99°E 


70.99°N 


296.68°E 


52.27°N 


348.24°E 


49.25°N 


329.42°E 



The time difference between the GEOSAT and SSM/I passes is less than 0.50 hours. 
The track of the GEOSAT through the SSM/I swath begins at 54.6°N,343°E and 
terminates at 65.4°N,324.4°E. The plot of SSM/I versus GEOSAT wind speeds (Fig. 
27) demonstrates a tremendous amount of scatter in the data. Wind speeds range from 
7.5 m/s to 20.5 m/s. The correlation calculated for this data set is -0.301. Looking at 
the plot of wind speed versus latitude (Fig. 28), one can see that there are two distinct 
regions where the difference in wind speeds are greater than 3 m/s. This occurs 
between 57.5°N and 60°N and again from 63.5°N until 65.3°N. The SSM/I wind field 
(Fig. 29) exhibits a range of wind speeds from 5 to 31 m/s. The highest wind speeds 
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occur along a band extending from the south central region of the field toward the 
coast of iceland. High wind speeds are also observed along the coast of greenland. 
The SSM/I water vapor field (Fig. 30) exhibits a high concentration of water vapor in 
the southeast quadrant of the SSM/I and GEOSAT orbits. The large white area in the 
water vapor image represents values higher than 60 kg/m 2 . An examination of the 
precipitation field (Fig. 31) discloses a major region of precipitation occurring along 
the track as well. The location of the precipitation coincides with the same area of 
Fig. 28 that reveals a large difference in wind speed between the SSM/I and GEOSAT 
altimeter (57.5°N to 60.0°N). The extent of precipitation appears to cover too small 
a region to be considered the only reason for the disparity between the GEOSAT and 
SSM/I wind speeds. The SSM/I cloud water field (Fig. 32) show's that this entire 
region is covered by high cloud water. From this evidence, the cloud water and 
precipitation together have apparently caused the large difference in wind speed values 
between the SSM/I and GEOSAT altimeter. 

A second area of concern lies between 63.5°N and 65.3°N (Fig. 28). This area 
is very close to the coast of Greenland, approximately 120 km at its nearest point. 
The altimeter wind speed measurement steadily increases as it approaches Greenland, 
while the SSM/I wind speed measurements decrease. This phenomena may be due to 
the closeness of the land mass and side lobe effects on the SSM/I. Criteria for 
nearness to land is not clear at this time. Chelton and McCabe [Ref. 1 8 j and Wentz 
et al [Ref. 19] consider only those regions where microwave remote sensing is done 
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SSM/ I VS ALTIMETER WIND SPEED 



ORBIT 7925 




ALTIMETER WIND SPEED C m/ 5 } 



Figure* 27 SSM/I vs. GEOSAT wind speed: orbit 7925. Line 1 
represents an actual correlation of -0.301, while line 2 represents a 
reference correlation of 1.0. 
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Figure 28 SSM/1 and GEOSAT wind speeds vs. latitude. 
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Figure 2 9 SSM/I wind speed field (m/s): orbit 7925. Line represents the 
GEOSAT flight path. 
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Figure 30 SSM/T water vapor field ikg/m 2 ): orbit 7925. Line through field 
represents the GEOSAT flight path. 
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Figure 31 SSM/1 precipitation field (mm/hr); orbit 7923. Line through field 
represents GEOSAT flight path. 
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Figure 32 SSM/T cloud water field (kg/m 2 ); orhit 7Q25. Line through field 
represents the GEOSAT flight path. 
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from 200 to 800 km from the nearest land mass. They sight erratic behavior due to 
side lobes and the change in emissivity. when making microwave measurements of the 
open ocean are made close to land. 
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VII. CONCLUSIONS 



To assess the relationship between the wind speeds as measured by the 
GEOSAT altimeter and the SSM/I, the correlation (R), R 2 , slope and intercept were 
calculated. Table one contains the results of the analysis of each data set. 

The value of R : represents the proportion of total variance in the values of the 
SSM/I wind speed data that can be accounted for by a linear relationship to the 
GEOSAT wind speeds. In orbits 9138 and 8058(6), more than 90% of the variance 
is explained by this linear relationship. Orbits 7925 and 8058(7) demonstrate almost 
no correlation between the SSM/I and altimeter wind speeds. 

The intercept gives the bias between wind speed measurements. Orbit 9138 
and 8058(6) exhibit the smallest bias, approximately 1 m/s. The accuracy goal of 



Table 1. Each data set is identified by its orbit number along with its 
associated correlation (R), R 2 , intercept, slope, and atmospheric parameters 
present during the wind speed measurement. 



ORBIT 


R 


R 2 


SLOPE 


INTERCEPT 


COMMENTS 


9138 


0.925 


0.8583 


1.284 


-1.297 


High water vapor 
concentration 


8934 


0.616 


0.5070 


0.644 


6.889 


Clear Atmosphere 


8058(7) 


0.208 


0.0431 


0.329 


10.599 


Precipitation, cloud water, and 
water vapor 


8058(6) 


0.935 


0.8739 


1.247 


-1.800 


Clear Atmosphere 


7925 


0.301 


0.0900 


0.258 


16.980 


Precipitation, cloud water, and 
water vapor 
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the GEOS AT altimeter and the SSM/I has been 2 m/s. Both orbits 9138 and 

8058(6) are within that goal. Orbits 7925, 8058(7), and 8934 have a bias in wind 
speed greater than 2 m/s. The slope is also a good indicator of the relationship 
between the SSM/I and GEOSAT altimeter wind speeds. A perfect relationship 
between wind speed measurements would give a slope of 1.0. Orbits 9138 and 
8058(6) have a slope very close to one, while 7925 and 8058(7) have a slope near 
0 . 0 . 

There are four conclusions that are made from the analysed data in this paper: 

(1) The data sets which did not contain cloud water or precipitation along the 

co-located SSM/1 and GEOSAT orbits (orbits 9138, 8934, and 8058(6)), 

demonstrated a good correlation between the two systems. 

(2) Although a large amount of water vapor was present along the entire track 
of orbit 9138, there appeared to be no evidence of bias with the measured wind 
speeds from either instrument. Since the SSM/1 corrects for the presence of water 
vapor in the atmosphere, it is obvious, from the high correlation between the 
GEOSAT altimeter and the SSM/I, that the SSM/I water vapor correction works 
well. The altimeter was not corrected for water vapor, suggesting that water vapor 
corrections on the altimeter were small. 

(3) The data sets which contained cloud water and precipitation (orbits 8058(7) 
and 7925), displayed large differences in recorded wind speeds between the SSM/I 
and GEOSAT altimeter. In these areas, the SSM/I wind speed would increase, 
while the altimeter wind speed decreased relative to one another and, consequently, 
very poor correlation resulted. There are several explanations for the difference in 
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wind speed relative to each other experienced by both sensors. The increased 
emissivity associated with precipitation, or cloud water, is the most likely cause of 
increased wind speeds experienced by the SSM/I, by introducing a "warm " 
atmospheric term. The GEOSAT may experience a decrease in wind speed at the 
same location. This may have been a result of a decrease in ocean surface 
roughness. The reduction in surface roughness is caused by rain drops impacting on 
the ocean surface and consequently attenuating the surface gravity waves [Ref. 21] . 
The reduction in surface roughness would result in a larger backscatter cross section 
for the GEOSAT altimeter. This reduced backscatter is interpreted as a decrease in 
wind speed, which is what is observed in the area of precipitation in orbits 7925 
and 8058(7). An alternative explanation for the low GEOSAT wind speed in the 
area of precipitation is that the measurement was correct and that the SSM/I had a 
bias to high wind speeds in the presence of rain. 

(4) Erratic behavior is seen in microwave measurements of the open ocean 
near a land mass. This is evident in orbit 7925 when measurements of wind speed 
and atmospheric properties were made within 120 km of the Greenland coastline. 
The resulting data from the coastal areas was highly erratic. Based on only one 
data set near a coastal feature, it appears that wind speed measured within 200 km 
of a land mass is not reliable. This is most likely due to the radiance detected by 
antenna sidelobes which mask the actual measurement of radiation from the main 
lobe of the antenna. 

The wind speed verification studies of the GEOSAT altimeter and the SSM/I 
demonstrate that both active and passive microwave systems measure similar values 
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of wind speed over the ocean. These instruments are comparable in open ocean 
regions at least 200 km from land and free of cloud water and precipitation. 
Utilization of GEOSAT altimeter and SSM/I wind speed data in a numerical 
forecasting model requires that a system for screening wind data for cloud water, 
precipitation, and proximity to land masses be used. 
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